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We show that a purely s-wave interaction in three dimensions (3D) can induce higher partial- 
wave resonances in mixed dimensions. We develop two-body scattering theories in all three cases of 
0D-3D, 1D-3D, and 2D-3D mixtures and determine the positions of higher partial-wave resonances 
in terms of the 3D s-wave scattering length assuming a harmonic confinement potential. We also 
compute the low-energy scattering parameters in the p-wave channel (scattering volume and effective 
momentum) that are necessary for the low-energy effective theory of the p-wave resonance. We point 
out that some of the resonances observed in the Florence group experiment [Phys. Rev. Lett. 104, 
153202 (2010)] can be interpreted as the p-wave resonances in the 2D-3D mixed dimensions. Our 
study paves the way for a variety of physics, such as Anderson localization of matter waves under 
p-wave resonant scatterers. 
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I. INTRODUCTION 

Scattering resonances play a central role in cold-atom 
physics [T]. In particular, s-wave Feshbach resonances 
induced by a magnetic field have been utilized to control 
the atom-atom interaction and have led to experimental 
realization of a rich variety of physics, such as the BCS- 
BEC crossover in Fermi gases [2 4 . Although a p-wave 
analog of the BCS-BEC crossover has been predicted the- 
oretically [SHS] and p-wave and higher partial-wave res- 
onances have been observed experimentally [U [TDHH], 
p-wave superfluids have not been realized in cold-atom 
experiments so far. This is because Fermi gases in the 
vicinity of the p-wave Feshbach resonances are unstable 
owing to inelastic collisions decaying into deeply bound 
dimers and do not reach their equilibrium within their 
short lifetime [20l [21]. In contrast, Fermi gases in the 
vicinity of the s-wave Feshbach resonances are long lived 
because such inelastic collisions are strongly suppressed 
by the Pauli exclusion principle (22j . 

In this paper, we propose a way to induce p-wave and 
higher partial-wave resonances from the s-wave Feshbach 
resonance only. This can be achieved by using a mix- 
ture of two atomic species A and B in mixed dimen- 
sions, where A atoms are confined in lower dimensions 
[quasi-zero dimension (OD), one dimension (ID), or two 
dimensions (2D)] and interact with B atoms living in 
three dimensions (3D) [231121] (see Fig. [l]). Indeed, the 
authors of Ref. [23 mention the possibility of p-wave res- 
onances in the 0D-3D mixture. Higher partial-wave res- 
onances induced from a purely s-wave interaction seem 
counterintuitive but can be understood by generalizing 
the argument for s-wave resonances given in Ref. [23L 

Suppose an s-wave interspecies interaction between A 
and B atoms supports a bound state with its binding 
energy Eb < 0. We confine only the A atom by a 3D 
harmonic potential with the oscillator frequency w. Be- 
cause the AB molecule also feels the confinement poten- 



tial, its center-of-mass motion is quantized. In the limit 
\Ei,\ Hoj or rriA ^ rriB, the energy of the AB molecule 
Eab is given by a sum of its center-of-mass energy and 
binding energy: 



Eab= {^+e + 2n 



rriA 



niA + mB 



■hLO-\Ebl (1) 



where mAirriB) is the mass of the A{B) atom, £ = 
0, 1, 2, . . . is an orbital angular momentum quantum 
number, and n > is an integer. We can view these 
quantized energy levels as a tower of £th partial-wave 
"Feshbach molecules." When one of them coincides with 
the A-B scattering threshold at ^huj, the coupling be- 
tween the center-of-mass and relative motions leads to 
the resonance occurring in the ^th partial-wave channel. 
Furthermore, there exists a series of resonances in each 
partial- wave channel corresponding to different values of 
n. 

Similarly, in the 1D-3D mixture, the A atom is confined 
by a 2D harmonic potential and the energy of the AB 
molecule is quantized as 



Eab = {l + \m\+2n) 



TTIA 



0D-3D 



mA + mB 



1D-3D 



■hLj-\Eb\, (2) 



2D-3D 



'Y? 



FIG. 1. (Color online) Illustrations of shallow p-wave 
molecules in mixed dimensions. Shown are the asymptotic 
angular wave functions of B atoms in 3D relative to A atoms 
confined in OD (origin), ID (z axis), and 2D [xy plane), re- 
spectively. See the text for details. 
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where m = 0,±l,±2,... is a magnetic quantum num- 
ber. When one of them coincides with the A-B scatter- 
ing threshold at hu, the resonance occurs in the jmlth 
partial-wave channel. There exists a series of resonances 
in each partial-wave channel corresponding to different 
values of n. 

In the 2D-3D mixture, the A atom is confined by a ID 
harmonic potential and the energy of the AB molecule is 
quantized as 



Eab - 



1 



niA 



■ huj — \Eh 



(3) 



When one of them coincides with the A-B scattering 
threshold at ^huj, the resonance occurs. Because the 
wave function of the AB molecule has the even (odd) 
parity for an even (odd) n, the associated resonance is 
in the even (odd)-parity channel, which is also phrased 
as an s-wave (p-wave) resonance in the 2D-3D mixture. 
There exists a series of resonances in each parity chan- 
nel corresponding to different values of n. In order to 
achieve these higher partial-wave resonances in mixed 
dimensions, one needs to apply a strong optical lattice 
only on A atoms without affecting B atoms [IS]. Such a 
species-selective optical lattice has been successfully im- 
plemented in a Bose-Bose mixture of ^^K and ®''Rb [26] 
and the 2D-3D mixed dimensions have been realized in 
the Florence group experiment [27\ 

The above simple arguments illustrate how the higher 
partial-wave resonances are induced in mixed dimensions 
from the purely s-wave interaction in a free space. In this 
paper, we develop two-body scattering theories in mixed 
dimensions, assuming the harmonic confinement poten- 
tial and the short-range interspecies interaction charac- 
terized by an s-wave scattering length (= 030) and ef- 
fective range (= tsd), which are extensions of analyses 
in Refs. [23l [24] aimed for s-wave scatterings only. The 
two-body scattering of a OD, ID, or 2D atom with a 3D 
atom is three-dimensional in the sense that it is described 
by three relative coordinates, and thus such scattering 
theories can be developed in parallel with the ordinary 
scattering theory in 3D. In particular, we determine the 
positions of all partial-wave resonances and will see that 
they are well understood from the above simple argu- 
ments in a wide range of the mass ratio uiA/mB ^ 1. 
We also compute the low-energy scattering parameters 
in the p-wave channel (scattering volume and effective 
momentum) that are necessary for the low-energy effec- 
tive theory of the p-wave resonance in mixed dimensions. 

Our analyses and results for the 0D-3D, 1D-3D, and 
2D-3D mixtures are presented in Sees. ]!!}]!!!] and |IV[ re- 
spectively, while some details are shown in Appendixes[A| 
[B] and[Cl respectively. Also in Sec jlV] we point out that 
some of the resonances observed in the Florence group 
experiment [37] can be interpreted as the p-wsve reso- 
nances in the 2D-3D mixed dimensions. Finally, Sec. |V| 
is devoted to a summary of this paper and the stabil- 
ity of confinement-induced molecules and their analogy 
with Kaluza-Klein modes in extra-dimension models are 



discussed. 

For later use, we define the reduced mass itiab = 
"m-Afng/ {mA + ti^), the total mass M = vtia + and 
the harmonic oscillator length /ho = \/hJrriAOJ. Below 
we set h = 1 and the range of integrations is assumed to 
be from —00 to 00 unless otherwise specified. 



II. 0D-3D MIXED DIMENSIONS 

A. Scattering theory [23J 

The scattering of a quasi-OD A atom with a B atom 
in 3D is described by a Schrodinger equation, 

2mA 2 ^ ^ 2ms y (4) 
= Eil){rA-,rB) 

for |ryi — r^l > 0. The short-range interspecies inter- 
action is implemented by the generalized Bethe-Peierls 
boundary condition [^ [55] : 



i^irA.rB) 



rATB^r 



1 1 

a{Ec) \rA-rB\ 



X(r). (5) 



Here a{Ec) is the energy- dependent "scattering length" 
defined by 



1 



~a(E, 



1 



- mABrsuEc 



(6) 



and the collision energy operator Ec in the present case 
is given by 



Er=E 



1 



2^2 



\ — mAijJ r 

2M 2 



(7) 



The solution to the Schrodinger equation Q can be 
written as 



i^{rA,rB) = MrA,rB) 
2n 



■niAB 



Jdr'GEirA, rB]r' , r')x{r') 



(8) 



where ipQ is a solution in the noninteracting limit and Ge 
is the retarded Green's function for the noninteracting 
Hamiltonian: 



GE{rA,rB]r'^,r'B) 



1 



{rA,rB\ 



E- Ho + ^:0^ 



\r'A,r'B) 



rriB 



rA 



(r'A) 



(9) 



e 



— \/2mB\J (nx-|-na-|-ni + |)w — B — iO+|rB— r^l 
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Here with quantum numbers n = (nxyUy^nz) is the 
normahzed wave function of an A atom in the 3D har- 
monic potential. 

We now consider the low-energy scattering in which 



2 2m B 

is satisfied, and then -00 becomes 

MrA,rB) = Ce'''-^^MrA), 



(10) 



(11) 



which represents the A atom in the ground state of the 
3D harmonic potential and the plane wave of B atom 
with the wave vector k. The asymptotic form of the 
wave function at a large distance S> Iho is given by 



ip(rA,rB) 



C 



+ 



-f{k,k') 



where f{k, k') with k' = kr-g defines the two-body scat- 
tering amplitude in the 0D-3D mixed dimensions: 



f{k,k') = - 



1 mB 
C mAB 



dr 



f ^ — ik' -r' 



0S(/)x(r'). (13) 



We note that x h&s an implicit k dependence. 

The unknown function % can be determined by substi- 
tuting the solution (|8| into the Bethe-Peierls boundary 
condition ([5|. Defining the regular part of the Green's 
function Q by 

GE(i'A,'t'B\r' ,r')\ 

— -6[r - r ) + Q{r;r ), 



we obtain 
1 



27r 
mAB 



dr'g{r;r')x{r'). 



(15) 

This integral equation determines x/C^ which in turn 
provides / from Eq. (13). 



Because of the 3D rotational symmetries of the system, 
/, X, and g can be decomposed into their partial-wave 
components: 



/(fc,fc')-E(2^+l)^^(^)^^(^-^')' 

i=0 
oo 

x(r)-^(2^+l)x^(r)P,(^-fc), 

oo 

g{r; r') = ^ (2^ + 1) g,(r; /) P,{f ■ f'). 



(16) 
(17) 
(18) 



Equations ( 13|) and lead to the tth partial- wave scat- 
tering ampl: 



itude given by 

1 TUB 



C niAB 



with 



~a{Ec) 



Xdr) = Cji{kr)(f>oir) 



27r 

rriAB 



(20) 



dr'ge{r;r')xe{r'). 



From the general argument based on the 3D rotational 
symmetries and unitarity [30j , or from the explicit calcu- 
lation that uses the Green's function in Eq. we can 
show that fi has the usual low-energy expansion: 



lim fe{k) = , (21) 

fe-o ' ' i^Wfc2 + o(fc4) + zfc2£+i' 



where a^g and r^g are effective scattering "length" and 
0o(''a), (12) "range" parameters in the fth partial-wave channel. Note 
that a|,g has the dimension of (length)^^+^ and r^^ has 



l-2£ 



(length) 

Substituting the expansion of fg into Eq. (19), we can 



determine the low-energy scattering parameters. In par- 
ticular, the effective scattering length a^g is given by 

(f) _ 1 TUB %/7r 



with 



CniAB 2^+1 



Xeir) = C 



dr'r"<i,l{r')xi{r') (22) 



'Mr) 



2^+1 (^+ i)! 

dr'gi{r;r')\i^^^Xi{r'). 



■rriAB 



(23) 



The £th partial-wave resonance in the GD-3D mixed di- 
mensions is defined by the divergence of a^^ — > oo, which 
occurs when 



1 



HE,) 
is satisfied. 



Xiir) 



27r 
rriAB 



dr'gi{r;r')\^^^Xe{r') (24) 



Positions of resonances 



We now solve the integral equation (24) numerically 



dr'j,ikr')<p*oir')Xiir') (19) 



to determine the positions of fth partial- wave resonances. 
For the purpose of illustrating qualitative results, we shall 
set r^Y) = 0. For quantitative predictions in a specific 
atomic mixture, it is necessary but straightforward to 
include the effective range correction [2Z1- Some details 
of our method to solve the integral equation are shown 
in Appendix [K\ 

Figure [2] shows the positions of the lowest five reso- 
nances in terms of Iho/o-sD for £ — 0, 1, 2, 3 partial- wave 
channels as functions of the mass ratio niA/rnB- For 
completeness, we have included the result for the s-wave 
0) resonance, which has been reported in Ref. |23j . 
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FIG. 2. (Color online) 0D-3D mixture: Positions of the lowest five resonances in terms of iho/iSD for I = Q (upper left), 1 = 1 
(upper right), 1 = 2 (lower left), and ^ = 3 (lower right) channels as functions of the mass ratio niA/mB- The dashed curves 
are from the approximate formula Eab = |fi<^ by using Eq. M with E^, = —fi^/{2mABaiTt). 



As we have discussed in Sec. |Tj there exists a series of 
resonances in each partial-wave channel induced from 
the purely s-wave interaction in a free space. Indeed, 
the resonance positions are well described by the ap- 
proximate formula Eab = by using Eq. ([l]) with 
Eb = —h?/{2mABa\D) hi a wide range of the mass ratio 
mA/mB ^ 1. For such a mass ratio, as is evident from 
Eq. ([T]), d-wave (/-wave) resonances are nearly degener- 
ate with s-wave (p-wave) resonances, and thus it could 
be difficult to distinguish them practically. We also note 
that the ^th partial- wave resonance is (2.^-1- l)-fold degen- 
erate in a spherically symmetric potential. This degener- 
acy is lifted when the confinement potential is deformed, 
and therefore the one resonance at a given Iho/o-su splits 
into more resonances. 



[Eq. ( 21 )] with keeping the two dominant terms at k 

1 



£ab 



ZniBa^g 



for e = o, 



for £> I. 



(25) 



Away from the resonance, these universal formulas are 
no longer valid. The binding energy eab = ~'«'^/(2to_b) 
has to be determined by solving the integral equation 



1 



a{Ec 



-xTir) 



27r 

rriAB 



dr'gi{r;r%^^^xTir'), (26) 



where xT ^ component of the spherical harmonics ex- 
pansion of X- 



X{r) 



oo z 



£^0 m^~. 



xT{r)Yr{f). 



(27) 



Confinement-induced molecules 



On the a^jj > side of every resonance, a shallow 
AB molecule is formed. In the vicinity of the resonance 



a;g > l^^ , its binding energy eab 



£; - fw < 



is determined by the pole of the scattering amplitude 



We note that the solution is independent of m, and there- 
fore there are 2^-1-1 degenerate molecules for a spherically 
symmetric potential. This degeneracy is lifted when the 
confinement potential is deformed. 

We now derive the asymptotic form of the molecular 
wave function in the vicinity of the ^th partial- wave reso- 
nance where /ho *C ^ |tb| is satisfied. From Eqs. (18]) 
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FIG. 3. (Color online) 0D-3D mixture: p-wave effective scattering volume ^off/^ho (upper figures) and effective momentum 
kculho (lower figures) as functions of iho/i3D for mass ratios vaA/mB ~ 6/40 (left), 1 (middle), and 40/6 (right). 



and ^ with the replacement k ^ in, we find 



1 niB e 



-ktb 



rr.=-l-- ■ -"'^^ (28) 

The angular parts of asymptotic wave functions of three 
degenerate shallow p-wave molecules ip ^ Y^'^^{rB) are 
illustrated in Fig. [T] 



D. Scattering parameters in the p-wave channel 

The effective scattering length in the s-wave channel 
acff = afg has been computed in Ref. [23]. Here we fo- 
cus on the p-wave (^ = 1) channel and determine its two 
low-energy scattering parameters, namely, the effective 



scattering volume VeB 



- ^(1) 

'off ■ 



and the effective momen- 



tum fcoff 

computed by eliminating C from Eqs. (|22| and (|23[) and 



The effective scattering volume can be 



solving the resulting integral equation numerically (see 
Appendix [a] for details) . In Fig. [s] Uoff /l^o for r^u — 
is plotted as a function of Iho/o-sD for three mass ratios 
itia/itlb = 6/40, 1, and 40/6. We confirm the existence 
of a series of p-wave resonances (vcs od) induced from 
the purely s-wave interaction in a free space, while they 
become narrower for larger ^ho/osD- We also find that 
the resonance is wider when a lighter atom is confined in 
lower dimensions. 



Similarly, the effective momentum can be computed 
from Eq. (19), and fceff^ho for ran = is plotted in Fig. [s] 



as a function of ^ho/isD for the same three mass ratios. 
In the vicinity of the p-wave resonance Vcs 3> Zhoi ''^eff 
and kcs determine the binding energy of three degener- 
ate shallow p-wave molecules via the universal formula 



(25). Both VeS and kes are important to the low-energy 



effective theory of the p-wave resonance discussed below. 



E. Low-energy effective theory 

The low-energy effective theory of the p-wave reso- 
nance in the 0D-3D mixed dimensions is provided by the 
action 



V2 \ 



(29) 



S = jdt^\{t) {idt)^A{t) 

-I- Jdt^lit) {idt + so)^j{t) 

+ goJdt [*^(t)V,*yi,0)<i>,(i) 

+$](i)V,*B(t,0)*^(t) 



where the summation over j — x,y, z is implicitly under- 
stood. 'i>A and '^B fields represent the A and B atoms in 
CD and 3D, respectively. The interaction between A and 
B atoms is described through their coupling with three 
p-wave molecular fields $j. go is their coupling strength 
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and So is the detuning from the resonance. These cut- 
off (A)-dependent bare parameters can be related to the 
effective scattering volume and effective momentum 
kcS by matching the two-body scattering amplitude from 
the action (29) with that shown in Eqs. (16) and (21). 



The standard diagrammatic calculation leads to the 
following scattering amplitude with collision energy s = 
fcV(2mB): 



iAik) 



ik ■ k' 



/{2mB)+eQ . niE 



By defining 



and 



9'o 



fA3 
K 3 



■ Afc2 + |ifc3) 



TUB 



2mBg^ Stt^ 



A EE 



■mB 1 



(30) 



(31) 



(32) 



we reproduce the scattering amplitude ( 21 ) in the p-wave 



1) channel up to a kinematical factor: 
27r 3fe • k' 



A{k) = -- 



rriB 



1 

■Ucff 



2 



fc2 



(33) 



This low-energy effective theory can be generalized easily 
to the case with more than one lattice site where B atoms 
are confined and could be used to investigate the many- 
body physics across the p-wave resonance as in the 3D 
case |31j . We note that an effective field theory of the p- 
wave resonance in 3D has been developed in connection 
with the a-n scattering [5^ I33j . 



F. Weak-coupling limit 

When asD < and |a3D| ^ ^ho, the confinement- 
induced resonances can be understood in a different way, 
as is discussed in Refs. [531 To the leading order in 
the weak-coupling expansion osD/^ho ~0, an A atom 
occupies the ground state in a 3D harmonic potential and 
creates a mean-field attractive potential felt by a i? atom. 
Therefore, the scattering of the B atom by the confined 
A atom is described by 



2m f 



27ra3D 
rriAB 



where \<t>o{r)\' = e'^' ''l^ /{^W,f . 



2ms 

(34) 

This Schrodinger 



equation, which is valid in the weak-coupling limit, is 
equivalent to the integral equation (15), where only the 
n = term is kept in the Green's function ^ and x(^) 
is identified as a-^Y:i(j)Q{r)il}{r). 



By matching the solution of Eq. ( 34 ) with the asymp- 



totic form (12), we can determine the scattering ampli- 



tude and low-energy scattering parameters in the weak- 
coupling limit IosdI ^ ^ho- In particular, the reso- 
nance occurs when a new bound state appears. This 



is possible even in the weak-coupling limit because 
the attractive potential becomes strong compared to 
the kinetic term by decreasing the mass ratio down 
to mA/niB ^ 1. We find that the resonances are 
achieved at the critical values of {ruB /mAB)ia3i)/lho) = 
-1.19,-7.89,-20.2,... in the ^ = channel [23l 
IMI, (ms/m^s) (flsDAho) = -5.36,-15.5,-31.2,... 
in the i ~ I channel, (mB /mAB){0'3D/lho) — 
-11.9,-25.6,-44.6,... in the £ = 2 channel, and 
(TOB/TOAs)(a3DAio) = -20.9,-37.9,-60.4,... in the 
i — i channel. 



III. 1D-3D MIXED DIMENSIONS 

A. Scattering theory [35j 

The scattering of a quasi-lD A atom with a B atom 
in 3D is described by a Schrodinger equation. 



V2 



^PA 1 2 2 ^PB 

-— + -m^w Pa - ^ - I ^gg^ 



X tP{pa;Pb,zab) = EiI;{pa,Pb,zab) 



for y^ipA ~ Pb) 



■"AB 



> 0, where p = (x, y), zab 



za — zb, and the center-of-mass motion in the z direction 
is eliminated. The short-range interspecies interaction is 
implemented by the generalized Bethe-Peierls boundary 
condition EHllMl: 



il{pA,PB,ZAB)\ 
1 



1 



^J{pA - PbY + z 



AB 



xip)- 



(36) 



The collision energy operator Ec in Eq. (|6| in the present 
case is given by 




2 2 

m^w p 



(37) 



The solution to the Schrodinger equation ( 35 ) can be 
written as 



i'iPA:PB,ZAB) = MPA,PB,Zab) 

2tt 



■niAB 



J dp'OsipA, Pb, zab; p' , p' ,0)xip'), 



(38) 



where ipo is a solution in the noninteracting limit and Ge 
is the retarded Green's function for the noninteracting 
Hamiltonian: 



GsiPA, Pb,zab; p'a^ p'bi^'ab) 

= (PA, PB:Zab\ 



'E- Ho + iO 

^TJlBmAB ^ , \ 
5- 2^(Pn{pA)4>n{PA) 



Pai PbiZab) 



(39) 



/2mB\/{nx+ny + l)i^-E-ia+\rB-r'g\ 



\rB 
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Here 4>n with quantum numbers n = {n^jUy) is the nor- 
mahzed wave function of an A atom in the 2D harmonic 
potential, and fs = {pb, — y''^^^ zab^ are coordinates 

of the B atom relative to the confined A atom. The 
anisotropic factor is such because separations in the x 
and y directions are associated with ms while a separa- 
tion in the z direction is associated with mAB [see the 
last two terms in Eq. (35 1]. 

We now consider the low-energy scattering in which 

(40) 



(41) 



E - uj = < uj 

2mB 

is satisfied, and then ipQ becomes 

i^oiPA,PB,ZAB) = Ce''"'''^0o(PA), 

which represents the A atom in the ground state of the 
2D harmonic potential and the plane wave of B atom 
with the wave vector fc. The asymptotic form of the 
wave function at a large distance {tbI S> Iho is given by 



iI'{Pa,Pb,zab) C 



rB 



-/(fc±,fei; 



(42) 

where f{k±,k'j_) with fc' = fcfs defines the two-body 
scattering amplitude in the 1D-3D mixed dimensions: 

1 



/(fei,fcl) = - 



C 



rriB 
mAB 



dp'e-^^'--P'^l{p')x{p') 



We note that x has an implicit fcj^ = [k^, ky) dependence 
and both x a-nd / depend on k through the Green's func- 
tion (pgl). 



The unknown function x can be determined by substi- 
tuting the solution ( 38 ) into the Bethe-Peierls boundary 



condition (36). Defining the regular part of the Green's 



function Q by 

Ge{pa, Pb,zab;p', p', 0) 

_ rriAB 



'^T^ViPA- Pb)^ + ^AB 

we obtain 

^ -X(p) = Ce^'=-^<^o(p) 



Pa,Pb^P,zai 

-S{P-P') 



Q{p;p') 



(44) 



27r 

mAB 



dp'g{p-p')x{p') 



(45) 

This integral equation determines x/C, which in turn 
provides / from Eq. ( 43 ) . 



Because the system has a rotational symmetry in the 
xy plane, /, x, and Q can be decomposed into their 
partial- wave components: 



Am arccos k± k , 



X(P)= E Xm(p)e™™''-^^ 

m— — oo 
oo 

g{p;p')= Sm{p;p')' 



^im arccos p- p 



(46) 
(47) 
(48) 



m— — oo 



Because m and —m are degenerate, we assume m > in 
this section without loss of generality. Equations ( 43 1 and 



( 45 1 lead to the mth partial- wave scattering amplitude 
1 



given by 
fm{k±, k'j_) = — - 



ms 
C V mAB 



with 



dp'j,n{k±p')(j)*o{p')Xni{p') 

(49) 



d{Ec) 



Xmip) = CJmik±p)(t)oip) 



271 



(50) 



dp'Qm{p;p')Xm.{p')- 



mAB 

From the explicit calculation that uses the Green's func- 



tion in Eq. 
low-energy expansion: 

lim fra{k) 



39 1 , we can show that has the following 



1 



0(k"l+'^k'Z\k"lk'j;'+'^) 



(51) 



where a"^^ and r^g-* are effective scattering "length" 
and "range" parameters in the mth partial-wave chan- 
(43) 

nel. Note that a^g-* has the dimension of (length) 2™+^ 



and r^g-* has (length) ^^2™. This unusual form of the low- 
energy expansion is owing to the lack of full 3D rotational 
symmetries. 



Substituting the expansion of into Eq. (49), we can 



determine the low-energy scattering parameters. In par- 
ticular, the effective scattering length a^^'' is given by 

1 



(2m -H)! („) 
(2™m!)2 



■niB 



C V mAB 2'"m! 



with 



dp'p"''<t>l{p')xM 
(52) 



Xni{p) = C 



1 



2™m! 
27r 

mAB 



p'^Mp) 



(53) 



dp'Gm{p-:P')\^^^Xni{p')- 



The mth partial-wave resonance in the 1D-3D mixed di- 
mensions is defined by the divergence of a"^'^ — )■ oo, 
which occurs when 



1 



~a{E, 
is satisfied 



Xm(p) 



27r 

mAB 



dp'g„(p;p')L^oX"('°') (54) 



Positions of resonances 



We now solve the integral equation ( 54 1 numerically to 



determine the positions of |m|th partial-wave resonances. 
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FIG. 4. (Color online) 1D-3D mixture: Positions of the lowest five resonances in terms of iho/asD for m — (upper left), m = 1 
(upper right), m = 2 (lower left), and m = 3 (lower right) channels as functions of the mass ratio niA/mB- The dashed curves 
are from the approximate formula Eab = fiu by using Eq. ^ with Ef, = —h^/{2mABa^o)- 



For the purpose of illustrating qualitative results, we shall 
set = 0. For quantitative predictions in a specific 
atomic mixture, it is necessary but straightforward to 
include the effective range correction [27]. Some details 
of our method to solve the integral equation are shown 
in Appendix [Bj 

Figure |4] shows the positions of the lowest five res- 
onances in terms of Iho/o-sD for m = 0,1,2,3 partial- 
wave channels as functions of the mass ratio mj^/mB. 
For completeness, we have included the result for the 
s-wave (m = 0) resonance, which has been reported in 
Ref. [24] . As we have discussed in Sec. |lj there exists a 
series of resonances in each partial- wave channel induced 
from the purely s-wave interaction in a free space. In- 
deed, the resonance positions are well described by the 
approximate formula Eab — by using Eq. ^ with 
£"{, = —h^/{2niAB0^Y)) in a wide range of the mass ra- 
tio TTLA/'mB ^ 1. For such a mass ratio, as is evident 
from Eq. ([2|, d-wave (/-wave) resonances are nearly de- 
generate with s-wave (p-wave) resonances, and thus it 
could be difficult to distinguish them practically. We also 
note that the |m|th partial- wave resonance for |m| > 1 
is twofold degenerate in a circularly symmetric potential. 
This degeneracy is lifted when the confinement potential 
is deformed, and therefore the one resonance at a given 
^ho/asD splits into two resonances. 



Confinement-induced molecules 



On the a"^^ > side of every resonance, a shallow 
AB molecule is formed. In the vicinity of the resonance 

a^g^ > 'hlr'"'^^' its binding energy eab = E - lo < 
is determined by the pole of the scattering amplitude 
[Eq. ( 51 )] with keeping the two dominant terms at A: — >■ 0: 



Sab 



0^ 



for 



0, 



1 



(55) 



(m) (m) 



for ItoI > 1. 



Away from the resonance, these universal formulas are 
no longer valid. The binding energy eab = — k^/(2tob) 
has to be determined by solving the integral equation 

1 27r f 

-Xmip) = / dp'Gmip; p')\^^^^Xmip')- (56) 



aiE, 



mAB 



We note that the solution for |m| > 1 is independent 
of the sign of m, and therefore there are two degener- 
ate molecules for a circularly symmetric potential. This 
degeneracy is lifted when the confinement potential is 
deformed. 

We now derive the asymptotic form of the molecular 
wave function in the vicinity of the jmlth partial- wave 
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FIG. 5. (Color online) 1D-3D mixture: p-wave effective scattering volume Vcs/l'to (upper figures) and effective momentum 
kculho (lower figures) as functions of iho/i3D for mass ratios vaA/rnB ~ 6/40 (left), 1 (middle), and 40/6 (right). 



resonance where Zho ^ ^ ^ l^^sl is satisfied. From 
Eqs. (38) and (39) with the replacement k — >■ in, we find 




4- i^-^y-prn^Q.I^B 



m— ±|t7i| £—\m\ 



2£+ 1 {l + m)\ 



mAB 



(57) 



One can see that different spherical harmonics ^ I7" 
with the same magnetic quantum number m contribute 
owing to the lack of full 3D rotational symmetries. The 
s-wave nature of the free-space interaction ensures that 
only £ = |to|, \m\ -|- 2, |77i| -I- 4, . . . contribute so that the 
wave function is an even function of zab- The asymp- 
totic behavior at a large separation Ifsj — > 00 is dom- 
inated by the i = \m\ component. The angular parts 
of asymptotic wave functions of two degenerate shallow 
p-wave molecules -0 ~ Y^^{rB) are illustrated in Fig. [l] 



computed by eliminating C from Eqs. ([52|) and (53) and 



solving the resulting integral equation numerically (see 
Appendix [b] for details). In Fig. [s] Woff/^ho fo'" ''3D — 
is plotted as a function of ^ho/osD for three mass ratios 
'm-A/'m-B = 6/40, 1, and 40/6. We confirm the existence 
of a series of p-wave resonances (ueS — >■ 00) induced from 
the purely s-wave interaction in a free space, while they 
become narrower for larger /ho/oao- We also find that 
the resonance is wider when a lighter atom is confined in 
lower dimensions. 

Similarly, the effective momentum can be computed 
from Eq. (49), and fccffZho for ''so = is plotted in Fig. [5] 



as a function of ^ho/<^3D for the same three mass ratios. 
In the vicinity of the p-wave resonance Vcg 3> /{Jq, foff 
and fcoff determine the binding energy of two degener- 
ate shallow p-wave molecules via the universal formula 
(55). Both Woff and k^s are important to the low-energy 



effective theory of the p-wave resonance discussed below. 



D. Scattering parameters in the p-wave channel 

The effective scattering length in the s-wave channel 
Ooff = a^g has been computed in Ref. [24 . Here we focus 
on the p-wave (|m| = 1) channel and determine its two 
low-energy scattering parameters, namely, the effective 
scattering volume VeB = Qeff ^^'^ effective momen- 
tum fcoff = r^Jj^. The effective scattering volume can be 



E. Low-energy effective theory 



The low-energy effective theory of the p-wave reso- 
nance in the 1D-3D mixed dimensions is provided by the 
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action 



J dtdz ¥^{t, z) (^dt + -^Ait, z) 



f ( V^ + V^' 

/ dtdzdp'fl(t,z,p) iidt + \_ " ) ■^B{t,z,p) 



2mB 



+ j dtdz z) (^dt + ^ + z) 

+ goJ dtdz [^lit, 2)V,*l5(t, z, 0)$,(i, z) 
+^]{t,z)Vj^Bit,z,0)^A{t,z) 



(58) 



where the summation over j = x,y is impUcitly under- 
stood. '^A and 'i>B fields represent the A and B atoms in 
ID and 3D, respectively. The interaction between A and 
B atoms is described through their coupling with two p- 
wave molecular fields (f>j. go is their coupling strength 
and Eq is the detuning from the resonance. These cut- 
off (A)-dependent bare parameters can be related to the 
effective scattering volume Vcs and effective momentum 
^eff by matching the two-body scattering amplitude from 
the action (58) with that shown in Eqs. (46) and (51). 



The standard diagrammatic calculation leads to the 
following scattering amplitude with collision energy e = 
P/{2mB): 



iA{k) 



ik I • fc'i 



fc^/(2mB)+eo I yrriBrnAE 



V 3 



Afc2 



By defining 



£0 ^ruBmAB _ ^JruBmAB 1 



9l 



37r2 



67r 



and 



y/mBmAB 
37r2 



A 



y/niBmAB fceff 

6Tr 2^' 



(59) 



(60) 



(61) 



we reproduce the scattering amplitude ( 51 ) in the p-wave 
(|m| = 1) channel up to a kinematical factor: 



Aik) 



2n 



3fe_L • fc'_L 



JmBrriAB — + ik^ ' 



(62) 



This low-energy effective theory can be generalized easily 
to the case with more than one tube where B atoms 
are confined and could be used to investigate the many- 
body physics across the p-wave resonance as in the 3D 
case [31]. We note that the low-energy effective theory 
of the s-wave resonance in the 1D-3D mixed dimensions 
has been derived and used to study three-body problems 
in Ref. m. 



F. Weak-coupling limit 

When asD < and {a^ul ^ ^ho, the confinement- 
induced resonances can be understood in a different way. 
To the leading order in the weak-coupling expansion 
aso/^ho — —0, an A atom occupies the ground state in 
a 2D harmonic potential and creates a mean-field attrac- 
tive potential felt by a i? atom. Therefore, the scattering 
of the B atom by the confined A atom is described by 



^ PB 



2m_B 
X iP{pb,zab) 



2mAB mAB 



\(I>o{pb)\'^5{zab) 



(63) 



2m_B 



1pipB,ZAB), 



where |0o(p)P = e-' ^ '■i" / {y/n l^o)'^ ■ This Schrodinger 
equation, which is valid in the weak-coupling limit, is 
equivalent to the integral equation (45), where only the 



n = term is kept in the Green's function (39) and x{p) 
is identified as a3D0o(p)V'(Pj 0)- 



By matching the solution of Eq. ( 63 ) with the asymp- 



totic form (42), we can determine the scattering ampli- 



tude and low-energy scattering parameters in the weak- 
coupling limit IcsdI ^ ^ho- In particular, the resonance 
occurs when a new bound state appears. This is possible 
even in the weak-coupling limit because the attractive 
potential becomes strong compared to the kinetic term 
by decreasing the mass ratio down to iriA/iTiB <C 1. We 
find that the reson ances are achieved at the critical val- 
ues of \/mB/mAB{a3i)/lho) = —0.730, —2.55, —4.34, .. . 

in the m — channel, \/mB/ rriAB (asD / ^ho ) — 
-1.96,-3.69,-5.44,... in the m — 1 channel, 
^J TTiB I rriAB (flsD/^ho) = -3.14,-4.85,-6.58,... in 
the m — 2 channel, and \/mB/ tuab (asD / ^ho) — 
—4.31, —6.00, —7.72, ... in the m — i channel. 



IV. 2D-3D MIXED DIMENSIONS 

A. Scattering theory | 35 | 

The scattering of a quasi-2D A atom with a B atom 
in 3D is described by a Schrodinger equation. 



^1 2 2 

-- — — + -rnyiaj Za — 

2myi 2 2niB 2m ab 



V2 



PAS 



(64) 



X ■4!{za,zb,Pab) EiI)(za,zb,Pab) 



for ^J{zA - zbY + p'ab > 0' where pab = {xa - 
xbtVa — Vb) and the center-of-mass motions in the x 
and y directions are eliminated. The short-range in- 
terspecies interaction is implemented by the generalized 
Bethe-Peierls boundary condition [28l [29] : 



iI^{za,zb.,Pab)\ 



1 



1 



j(^c) \/{zA - zbY + Pab 



X{z)- 



(65) 
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The collision energy operator Ec in Eq. (|6|) in the present 
case is given by 



V2 



1 



2^2 



2M + 2"''-^" ' 



(66) 



The solution to the Schrodinger equation (64 1 can be 
written as 



tp{zA,ZB,PAB) = iPo{za,zb,Pab) 

+ 

VtlAB 



dzCsizA, zb^Pab; z', z', 0)x{z'), 



(67) 



where ipQ is a solution in the noninteracting limit and Ge 
is the retarded Green's function for the noninteracting 
Hamiltonian: 



Ge{za, zb,Pab', z'a, z'btP'ab) 
1 



{za,zb, Pab\ 



E-Ho + iO^ 



I^'aj z'btP'ab) 



(l^nA^AWnA^A) 



mAB sr^ 
' 27r 

^-^2mB\/{n^ + ^)u:-E-iO+\fB-f'i^\ 



(68) 



\rB-r'g\ 

Here </)„^ is the normalized wave function of an 
A atom in the ID harmonic potential, and rs = 

^ ^-"^^ Pabi zb^ are coordinates of the B atom rel- 
ative to the confined A atom. The anisotropic factor is 
such because a separation in the z direction is associated 
with niB while separations in the x and y directions are 
associated with mAB [see the last two terms in Eq. (64)]. 
We now consider the low-energy scattering in which 



E--UJ. 
2 



fc2 



2tob 

is satisfied, and then -00 becomes 

iI^o{za,zb,Pab) ^ Ce^^"^''(j}a{zA), 



(69) 



(70) 



which represents the A atom in the ground state of the 
ID harmonic potential and the plane wave of B atom 
with the wave vector k. The asymptotic form of the 
wave function at a large distance Ifsj 3> ^ho is given by 



^{za,zb,Pab) ^ C e*'=-^« + -: — f{kz,K) M^a), 

I I^B J 

(71) 

where /(fc^, k'^) with k' = kf-B defines the two-body scat- 
tering amplitude in the 2D-3D mixed dimensions: 



1 

C 



dz'e 



0S(z')x(^')- (72) 



We note that x ti&s an implicit k^ dependence and both 
X and / depend on k through the Green's function (68 1. 



The unknown function x can be determined by substi- 
tuting the solution (67) into the Bethe-Peierls boundary 



condition ( 65 ) . Defining the regular part of the Green's 



function Q by 

Ge{za,zb,Pab;z',z\0) 



rriAB 



za,zb^z;Pab^O 

6iz-z') + g{z;z'), ^^^^ 



2t^\/{za - zbY + Pab 
we obtain 

^ X{z) = Ce^'^'^Uz) + — / dz'g{z- z')x{z'). 



~a{Ec) 



mAB 



(74) 

This integral equation determines x/C*, which in turn 
provides / from Eq. ( 72 ) . 



Because the system has a refiection symmetry about 
the z axis, /, X: and G can be decomposed into their 
even- and odd-parity components: 

/±(fc.,fc^)= ^^^^'^^^y^^^'~^^\ (75) 



G±{z;z') 



X{z) ± x{-z) 
2 

giz;z')±giz;~/) 



(76) 
(77) 



Equations ( 72 ) and ( 74 ) lead to the even- and odd-parity 



scattering amplitudes given by 

-ik^z' _j_ f^k^z 



f± [kz, k^ 
with 



(78) 



1 pik^z I -ik^z 

x±{z)^c Mz) 



d{E,) 



+ 



2n 

mAB 



(79) 



dz'g±{z:z')x±{z'). 



From the explicit calculation that uses the Green's func- 



tion in Eq. ( 68 1 , we can show that f± has the following 



low-energy expansion: 
lim /+(fc^,A:;) 

A;— 



_i i^(+) 

,(+) 2 

^off 



r^/P o{k^) +ik[l + Cl(fc2)] 



(80) 



0{klk'^' 



and 



lim /_(A;^,fc;) 



3 



^ - ir^ff'fc2 + o(fc4) + ifc3 [1 + 0{k2)] 



(81) 



+ o{ex,kzkf), 



where and r^g"* are effective scattering "length" and 
"range" parameters in the even- and odd-parity chan- 
nels. Note that has the dimension of (length)^^-'^ 
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FIG. 6. (Color online) 2D-3D mixture: Positions of the lowest five resonances in terms of /ho/asD for even-parity (left) and 
odd-parity (right) channels as functions of the mass ratio ttia/ttis. The dashed curves are from the approximate formula 
Eab = by using Eq. (jsl with Ei, = -h^ /{2mABalj^). 



and r^g-* has (length)*^. This unusual form of the low- 
energy expansion is owing to the lack of full 3D rotational 
symmetries. 



Substituting the expansion of f± into Eq. ( 78 ) , we can 



determine the low-energy scattering parameters. In par- 
ticular, the effective scattering length is given by 



1 



~a{E,) 



X+{z) = -^(l)o{z) I dz' (f)l{z')x+{z') 



+ 



2% 
mAB 



(82) 



and 



X-(^) ^ — 7^z(t)o{z) / dz'z'<j)'^{z')x-{z') 



3a 



(-) 
27r 
mAB 



(83) 



dz'g^{z;z')\i^^^X-{z'), 



where we have eliminated C. The even- and odd-parity 
resonances in the 2D-3D mixed dimensions are defined 
by the divergence of a|,^'' — )■ oo, which occurs when 

r^X±W = ^Jdz'gAz;z%^,x^{z') (84) 

is satisfied. 



Figure [6] shows the positions of the lowest five res- 
onances in terms of /ho/isD for even- and odd-parity 
channels as functions of the mass ratio toa/^^b- For 
completeness, we have included the result for the even- 
parity resonance, which has been reported in Ref. |24) . 
As we have discussed in Sec. [l] there exists a series of res- 
onances in each parity channel induced from the purely 
s-wave interaction in a free space. Indeed, the resonance 
positions are well described by the approximate formula 
Eab = by using Eq. ^ with Et = -h^ / {2mABai^) 
in a wide range of the mass ratio mAlvriB ^ 1- 



C. Confinement-induced molecules 

On the a^g-* > side of every resonance, a shallow 
AE molecule is formed. In the vicinity of the resonance 



> its binding energy eab 



E 



< 



is determined by the pole of the scattering amplitude 
[Eqs. (80 1 and (81 )] with keeping the two dominant terms 



at fc -> 0: 



Sab 



2mBa 
1 



(+)2 
cff 



for even parity. 



for odd parity. 



(85) 



Positions of resonances 



We now solve the integral equation (84) numerically 



to determine the positions of even- and odd-parity reso- 
nances. For the purpose of illustrating qualitative results, 
we shall set — 0. For quantitative predictions in a 
specific atomic mixture, it is necessary but straightfor- 
ward to include the effective range correction [27]. Some 
details of our method to solve the integral equation are 
shown in Appendix [C| 



Away from the resonance, these universal formulas are 
no longer valid. The binding energy eab = / {2mB) 
has to be determined by solving the integral equation 

-^,X±iz) = — ldz'g^{z-z')l^^^X±{z'). (86) 

We now derive the asymptotic form of the molecular 
wave function in the vicinity of the even- or odd-parity 
resonance where Zho ^ ^ ItbI is satisfied. From 
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FIG. 7. (Color online) 2D-3D mixture: p-wave effective scattering volume Woff/^ho (upper figures) and effective momentum 
fccft^ho (lower figures) as functions of iho/i3D for mass ratios vaA/mB ~ 6/40 (left), 1 (middle), and 40/6 (right). 



Eqs. (67) and (68) with the replacement k — >■ zk, wc find 



I— even or odd B 



xY^ifs) dz'\z'\'<t>l{z')x±{z'). 



(87) 



One can see that different spherical harmonics ^ con- 
tribute owing to the lack of full 3D rotational symme- 
tries. The s-wave nature of the free-space interaction en- 
sures that only the m = component contributes so that 
the wave function is independent of Pab- The asymp- 
totic behavior at a large separation jr^l oo is domi- 
nated by the ^ = or 1 component in the even- or odd- 
parity channel, respectively. Therefore, we can phrase 
the even (odd)-parity resonance as an s-wave (p-wave) 
resonance in the 2D-3D mixture. The angular part of the 
asymptotic wave function of a shallow p-wave molecule 
■0 ~ y]°(fB) is illustrated in Fig. [l] 



scattering volume can be computed by solving the inte- 
gral equation J83[ ) numerically (see Appendix [C] for de- 
tails). In Fig. ^n;off/^hQ for rso = is plotted as a func- 
tion of ^ho/flsD for three mass ratios mA/mB = 6/40, 1, 
and 40/6. We confirm the existence of a series of p-wave 
resonances {vcs — oo) induced from the purely s-wave 
interaction in a free space, while they become narrower 
for larger /ho/asD- We also find that the resonance is 
wider when a lighter atom is confined in lower dimen- 
sions. 

Similarly, the effective momentum can be computed 
from Eq. (78), and fccffZho for ran = is plotted in Fig. [t] 



as a function of ^ho/<^3D for the same three mass ratios. 
In the vicinity of the p-wave resonance VcS 3> /jj^, Vcs 
and fceff determine the binding energy of a shallow p- 
wave molecule via the universal formula (85). Both VeS 



and fcoff are important to the low-energy effective theory 
of the p-wave resonance discussed below. 



D. Scattering parameters in the p-wave channel 

The effective scattering length in the s-wave (even- 
parity) channel acs = a^g'' has been computed in 
Ref. ;24] . Here we focus on the p-wave (odd-parity) chan- 
nel and determine its two low-energy scattering parame- 
ters, namely, the effective scattering volume VeS = a^g'' 
and the effective momentum fccff = '"Iff^- The effective 



E. Low-energy effective theory 



The low-energy effective theory of the p-wave reso- 
nance in the 2D-3D mixed dimensions is provided by the 
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action 



F. Comparison to experiment | 27| 



+ / dtdpdz '^\{t,p,z)\idt+ \-^B{t,p,z) 

+ j dtdp<^\t,p) ^idt 
+ go dtdp 



vi/^(i,p)v,*yi,p,o)$(t,p) 

$^(i,p)V,4's(i,p,0)*A(t,p) 



^ffyi and ^ffi fields represent the A and B atoms in 2D 
and 3D, respectively. The interaction between A and B 
atoms is described through their coupling with a p-wave 
molecular field $. go is their coupling strength and Eq 
is the detuning from the resonance. These cutoff (A)- 
dependent bare parameters can be related to the effec- 
tive scattering volume VcS and effective momentum fceff 
by matching the two-body scattering amplitude from the 
action (88) with that shown in Eqs. (751 and (81). 



The standard diagrammatic calculation leads to the 
following scattering amplitude with collision energy e — 
e/{2mB): 



iA{k) 



By defining 



and 



1 kz k 



9^0 



37r2 



(A^ 
\ 3 



£0 ^ ruAB ^ rriAB 1 



i l_ mAB 1^ ^ ruAB fccff 



2mBgQ 37r-^ 



(89) 



(90) 



(91) 



we reproduce the scattering amplitude ( 81 ) in the p-wave 



(odd-parity) channel up to a kinematical factor: 



A{k) = - 



27r 



3/t y A^^ 



mAB — + ik^ ' 



(92) 



This low-energy effective theory can be generalized easily 
to the case with more than one layer where B atoms 
are confined and could be used to investigate the many- 
body physics across the p-wave resonance as in the 3D 
case [31]. We note that the low-energy effective theory 
of the s-wave resonance in the 2D-3D mixed dimensions 
has been derived and used to study many-body problems 
in Refs. [371 [3H]. 



Finally, we compare our predictions to the experi- 
mental measurement of resonance positions reported in 
Ref. [13 . The Florence group has realized the 2D- 
3D mixed dimensions by using a Bose-Bose mixture of 
A ='*^K and B =^^Rb with a species-selective ID optical 
lattice: 



^2 

VizA) = s-^sm'^{kLZA). 
ZniA 



(93) 



Here s is the lattice depth parameter and k^ = 2tt/Xl 
with Al = 790.02 nm is the wave vector of the laser 
light. By monitoring three-body inelastic losses, a series 
of resonances has been observed as a function of the mag- 
netic field. The measured resonance positions are shown 
in Fig. [8] for five values of the lattice depth parameter, 
s = 10, 15, 20, 23, and 25. 

As for theoretical predictions, we assume that the lat- 
tice potential ( 93 ) can be approximated by a harmonic 
potential: 



V{za) 



1 



(94) 



The harmonic oscillator length is given from the laser 
wavelength by /ho = AL/(27rs^/*). With the use of the 
free-space value of the effective range r^u = 168.37 ao 
for the ■^^K-^^Rb mixture, we solve the integral equation 



( 84 ) numerically to determine the resonance positions in 
terms of /ho/osD- The free-space scattering length 
thus obtained is converted into the magnetic field B (G) 
by the following empirical formula |27| : 



ao 



= 208 1 



30.9 



49.92 



1.64 



B + 38.52 B- 38.37 S - 78.67, 

(95) 

where ao is the Bohr radius and the resonance in a free 
space occurs at B — 38.4 G. These critical magnetic field 
values for the s-wave (aos — >■ oo) and p-wave {vcS oo) 
resonances are plotted in Fig. [8] as functions of the lattice 
depth parameter s. 

In Fig. [8) the resonances are labeled by an integer 
n — 0,1,2,3 in descending order, corresponding to n 
m Eq. ^. One can see the reasonable agreement be- 
tween the experimental measurements and the theoreti- 
cal predictions for n — (s-wave) and n = 1 (p-wave) 
resonances. They begin to deviate for larger n simply 
because the harmonic potential approximation ( 94 ) be- 



comes worse for higher excited states in an optical lattice. 
Although the lattice potential (93) has to be taken into 
account to improve the quantitative agreement [27], for 
a sufficiently strong optical lattice, the resonances cor- 
responding to odd n should be the p-wave resonances 
in the 2D-3D mixed dimensions, as we have discussed 
in this section. Further experimental investigations to 
confirm the p-wave nature of these resonances would be 
worthwhile. The p-wave nature of shallow molecules as 
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10 15 20 25 30 

lattice depth parameter s 

FIG. 8. (Color online) Experimentally measured resonance 
positions in terms of the magnetic field B as functions of the 
depth parameter s of an optical lattice [27j. Also shown are 
the predicted positions of s-wave (aoff — > oo; solid curves) and 
p-wave (uofi — > oo; dashed curves) resonances in a harmonic 
potential as well as the resonance position in a free space 
(osD — >■ 00 ; dotted line). 

illustrated in Fig. [T] could be seen as a suppression of the 
density of B atoms, as opposed to an enhancement for 
s-wave molecules, at the center of the confinement po- 
tential after a sudden change of v^s > to its negative 
side. 



V. SUMMARY AND CONCLUDING REMARKS 



In this paper, we showed that a purely s-wave interac- 
tion in a free space can induce higher partial-wave res- 
onances in mixed dimensions. We developed two-body 
scattering theories in all three cases of 0D-3D, 1D-3D, 
and 2D-3D mixtures, and determined the positions of 
higher partial- wave resonances in terms of the free-space 
scattering length, assuming a harmonic confinement po- 
tential. We also computed the low-energy scattering pa- 
rameters in the p-wave channel (effective scattering vol- 
ume and momentum) that are necessary for the low- 
energy effective theory of the p-wave resonance. Poten- 
tially our study paves the way for a variety of physics, 
such as Anderson localization of matter waves under 
p-wave resonant scatterers in the 0D-3D mixed dimen- 
sions [231 ESI ED], Shiba bound states in a superfiuid 
medium |34) , and Bose-Einstein condensation of p-wave 
molecules in the 2D-3D mixed dimensions, provided that 
such molecules are long lived. 

When both A and B atoms are fermionic (such as for 
the ^Li-''^°K mixture [4TlI45j ) . inelastic three-body, atom- 
molecule, and molecule-molecule collisions decaying into 
deeply bound dimers whose size is set by the range of 
interatomic potential tq are suppressed. This is because 



at a short distance '-^ ro <C /ho, the confinement poten- 
tial is irrelevant, and in the s-wave interspecies interac- 
tion the Pauli exclusion principle is effective to suppress 
the inelastic collisions decaying into the deeply bound 
dimers (52} . An order-of-magnitude estimate of such a 
three-body recombination rate can be found in Refs. |36j 
and for wide Feshbach resonances and in Ref. [55] for 
narrow Feshbach resonances. 

Therefore, for a Fermi-Fermi mixture in mixed dimen- 
sions, the relaxation of molecules is dominated by the 
decay into deeper molecular states of size ^ /ho if they 
exist. In the vicinity of the s-wave resonance at the small- 
est value of /ho/asD, indicated by thick curves in Figs. [2] 
|4j and[6j there is no such molecular state and thus asso- 
ciated s-wave molecules are long lived. However, in the 
vicinity of the higher partial-wave resonances, there is 
always at least one deeper molecular state. The inelas- 
tic atom-molecule and molecule-molecule collisions de- 
caying into such molecular states are not generally sup- 
pressed, and thus p-wave or higher partial- wave molecules 
in mixed dimensions would be short lived. It is therefore 
an important future problem to investigate as to whether 
the inelastic collisions of p-wave molecules decaying into 
the lowest s-wave state of size /ho could be suppressed, 
for example, by controlling the system parameters. 

We conclude this paper by discussing a simple analogy 
between a tower of our confinement-induced molecules 
and Kaluza-Klein modes in extra-dimension models [IB] . 
Suppose our 3D world has one compact extra dimension 
with an extent of L. Then the momentum of a 4D particle 
in the extra fourth direction is quantized as p4 = i-nhn/ L. 
Because we live in 3D, such a particle can be viewed as 
a tower of "new particles" with the quantized masses 
TO„ — 2TThn/cL. These new particles are called Kaluza- 
Klein modes and could be observed by colliding high- 
energy particles > 2TThc/L at the Large Hadron Collider. 

Now in our case, for example, in the 2D-3D mixed 
dimensions, the motion of a 3D molecule in the "ex- 
tra" z direction is quantized as in Eq. ^ because of 
the confinement potential. We can view such quantized 
energy levels as a tower of new "Feshbach molecules" 
analogous to the Kaluza-Klein modes. Unlike in extra- 
dimension models, our particle is a composite molecule 
and we can control its binding energy to shift the ener- 
gies of the tower of new molecules. When each of them 
crosses the scattering threshold, the resonance occurs. A 
series of such resonances caused by the "Kaluza-Klein" 
tower of confinement-induced molecules has been inves- 
tigated in this paper and observed in the cold-atom ex- 
periment [531157) . 
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Appendix A: 0D-3D mixed dimensions 



Here we present some details of our method to solve integral equations derived in Sec. |ll] The central issue is the 

3, , _ 
2^ = 2WrE 



evaluation of the regular part of the Green's function defined in Eq. |l4| . When E — '^uj = < 0, it is useful to 
represent Q as 

g(r;r') 



dTe^"^B — ( exp 

\27rsinhwT/ \2nT/ 



niAi^ [r + r') coshwr — 2r • r' 



l\2 



sinhc 



2t 



(r - r') 



(Al) 



Its partial-wave projection gives 



gi{r;r') = -J dcos9g{r;r') PAcosi 



X exp 



= / m^o^e-- / 7715x3/2 . 

dre^'-B I- — — - — i ji 

\ ZiT smh UJT J \ Zttt / 

f mAio coshiOT TUB 



mAu: mB\ I 

I I 1 rr 

smhwr r / 



V sinh UJT 



dr 



^AB \ '^/^ 5{r — r') 



(A2) 



We now evaluate the matrix elements of Qg with respect to the eigenfunctions of 3D harmonic oscillator with orbital 
angular momentum t. 



1 



;3/2y2.(„ + ^+i). 



'ho 



(A3) 



which form an orthonormal basis: 



dr 0W {r)^f} (r) = / Anr^dr 0(f) (r)0i? (r) = 6„ 



(A4) 



A lengthy but straightforward calculation leads to 



Mff ^ ^ / drdr'cj,f\r)ge{r;r')^f\r') 



rriAB 



-1 Fmi TUB fi + l+h\fj + l + h\ f°° ^dr f e"^ ^ 



40r V mA niAB 



X 1 



X 



A2 -r2 



2F1 



j J Jq (a;T)3/2 ysinhwT 

2" 



1 



Y 



2' V A2 - A-r2 



mAB 



A2-r2 

Lodr 



2TrmA Jo (a;T)3/2 • 



(A5) 



where 2F1 is the hypergeometric function and 



X = 



1 / 1 



TTlB 



2 V tanhwr tm^wt 



and Y = 



If 1 



777b 



2 V sinhwr niAUOT 



(A6) 



Expanding xi in terms of 4>n\ the integral equations in Sec. |ll] reduce to linear algebraic equations that have been 
solved numerically. 
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Appendix B: 1D-3D mixed dimensions 

Here we present some details of our method to solve integral equati ons derived in Sec. |III[ The central issue is the 
valuation of tl 
represent Q as 



evaluation of the regular part of the Green's function defined in Eq. (|44|). When E — uj = < 0, it is useful to 



f J mAOJe'^^ niB niAB 

^'^''^'^ io 27rsinha;T 27rT V Sttt 

Its partial-wave projection gives 

Gm{p;p') = - difQ{p\p')cos{mLp) 

Jo 



exp 



m^o; (p + p')^ coshwT 2p ■ p' ms 



sinh L 



2t 



[P-P'f 



(Bl) 



X exp 



, -h^r mAUJe^' rriB mAB ^ 
dre^^B — ./ 7„ 

27rsinhaJT 2'kt y 2ttt 

/ mA(^ COShujT TTis^ + p'^ 

\ sinh UJT T 



Y 2ttt ™ L \ sinh ujt t J 



UJT 

^^^(mABY^^ Sjp-p') 
^ V 27rr / 27rp 



(B2) 



We now evaluate the matrix elements of Q^^ with respect to the eigenfunctions of 2D harmonic oscillator with 
magnetic quantum number m: 



4rHp)^^J- 



(n + m)! 



ho 



LT^p'/iL), 



which form an orthonormal basis: 



dp '^i"^ (P).^!"^ (p) = 27rp dp (p) = 5^ 

A lengthy but straightforward calculation leads to 



(B3) 



(B4) 



(m) _ 27r/ho 



mAB 
-1 niB 

2y/2TT ^JmAmAB 
X 1 



dpdp'4"^\p)g„.ip-p')<l>'r\p') 



i + m\ fj + m 



X2 -y2 



J /Jo 
i,-j;m+ 1 



(a;r)3/2 



Y 



X^-X-Y^ 



sinhwT 

2" 



X2 -y2 



Lodr 



(B5) 



27rm^ 7q (wt)3/2 



Expanding Xm in terms of ^n"'', the integral equations in Sec. Ill reduce to linear algebraic equations that have been 
solved numerically. 



Appendix C: 2D-3D mixed dimensions 



Here we present some details of our method to solve integral equations derived in Sec. |IV[ The central issue is the 



evaluation of the regular part of the Green's function defined in Eq. (73 1. When E - 
represent Q as 

Q{z;z') = - 



< 0, it is useful to 



dr t 



rriB mAB 



27rsinhtjr V 2nT 2ttt 



exp 



TOylW (z + z')2 cosh WT — 2zz' TTlB 



sinhwr 



2r 



(Cl) 
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Its parity projection gives 



dr ( 



X exp 



TOAwe""^ ImB niAB fcoshl 
27r sinh V 27rr 2-KT 1 sinh J 

/ myii^ cosh cjr z'^ + z' 

\ sinh wr t 



^^i^ ^ I 

smh T 

mAB\^''^ 5{z~z')±5{z + z') 
2ttt 



We now evaluate the matrix elements of Q± with respect to the eigenfunctions of ID harmonic oscillator: 

1 ^-1/4 



(jjniz) 



\/2"n! 



which form an orthonormal basis: 



dz(j)l_{z)(j)n-{z) = Sn„>- 



A lengthy but straightforward calculation leads to 

^(+) _ 27r/ho Uzdz'^,{z)g+{z-z')cpj{z') 



2V'rV ™A (f )!(!)! V 7o (wt)3/ 



Lodr 



X 1 



2' 2'2'VX2-X-r2 



sinh WT 

2" 



X2_y2 



1/2 



LodT 



2-KmA Jo (wt)3/2 



for i,j = even integers and 



(-) 



2ttI 



ho 



dzdz' (j)i{z)Q^{z; z')(j)j{z' 
1 /m^ (—1)3 + 5 / i! j\ 



^V"^a(^)!(^)!V2'+j7o (^t)3/ 



ujdr 



1 



X 1 



X2 -y2 



2F1 



I- 1 i-j 3 



sinhwr ^ - 

2^ 



3/2 



Y 



2 '2'Vx2-x-y2 



2'iTmA Jo (wr)3/2 



(C2) 



(C3) 



(C4) 



(C5) 



(C6) 



for i,j — odd integers. Expanding x± i^i terms of 0„, the integral equations in Sec. IV reduce to linear algebraic 
equations that have been solved numerically. 
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